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Abstract 

We develop an extension of Bohmian mechanics by defining Bohm-fike tra- 
jectories for (one or more) quantum particles in a curved background space-time 
containing a singularity. Part one, the present paper, focuses on timelike singu- 
larities, part two will be devoted to spacelike singularities. We use the timelike 
singularity of the (super-critical) Reissner-Nordstrom geometry as an example. 
While one could impose boundary conditions at the singularity that would pre- 
vent the particles from falling into the singularity, in the case we are interested in 
here particles have positive probability to hit the singularity and get annihilated. 
The wish for reversibility, equivariance and the Markov property then dictate that 
particles must also be created by the singularity, and indeed dictate the rate at 
which this must occur. That is, a stochastic law prescribes what comes out of the 
singularity. We specify explicit model equations, involving a boundary condition 
on the wave function at the singularity, which is applicable also to other versions 
of quantum theory besides Bohmian mechanics. 

Key words: quantum theory in curved background space-time; Reissner-Nordstrom 
space-time geometry; timelike singularities; Bohmian trajectories; particle cre- 
ation and annihilation; stochastic jump process. 

1 Introduction 



We study the construction and behavior of quantum mechanics in a space-time with 
(timelike) singularities, using Bohmian mechanics (also known as pilot-wave theory), a 
version of quantum mechanics with particle world lines. We develop an explicit model, 
defining trajectories on a given background space-time geometry, using the Dirac equa- 
tion, Bohm's law of motion, and a novel boundary condition at the singularity that 
governs the creation and annihilation of particles at the singularity. As a case study, 
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we consider the Reissner-Nordstrom geometry in the charge > mass regime, featuring 
a naked timehke singularity. 

Bohmian mechanics was developed as a realist version of nonrelativistic quantum 
mechanics [1] and succeeds in explaining all phenomena of quantum mechanics in terms 
of an objective, observer-independent reality consisting of point particles moving in 
space; see [13] for an overview. Bohmian mechanics possesses a natural extension to 
relativistic space-time if a preferred foliation of space-time into spacelike hypersurfaces 
(called the time foliation) is granted 0. This extension has also been formulated for 
curved space-time geometries [25l [26] , but not yet for geometries with singularities. 
While horizons present no difficulty, singularities require further work to define the 
theory: Basically, we have to specify what happens when a particle hits the singularity, 
since at this point the law of motion is no longer defined. The possibility we consider 
here is that the particle gets annihilated: that is, if the system consisted of N particles, 
then the world line of the particle that has arrived at the singularity ends there, while 
the other — 1 particles, which are not at the singularity and thus have no reason to 
vanish, continue to move according to Bohm's law of motion. To make this possible, 
we need wave functions from Fock space, as always when particles can get created or 
annihilated. Further considerations then naturally lead us to specific equations, defining 
a Bohm-type theory. 

1.1 Particle Creation 

The most astonishing feature of this theory is perhaps that particles can also be created 
at the singularity. Indeed, this occurs at random times distributed in a way determined 
by a law involving the wave function. As a consequence, the theory is no longer de- 
terministic; instead, the evolution of the particles can be described by a Markovian 
stochastic process in the appropriate configuration space. This process has similarities 
with processes used earlier in "Bell-type quantum field theories" [TU] , a stochastic version 
of Bohmian mechanics for quantum field theories that we describe briefly in Section [61 

It is an interesting moral from our model that particle creation and annihilation, 
one of the characteristic features of quantum field theory, can arise from the presence 
of timelike singularities. This invites speculation whether quantum particles capable of 
emitting other particles might actually contain timelike singularities. After all, the clas- 
sical gravitational field of a point particle with the charge and the mass of an electron is 
exactly the Reissner-Nordstrom geometry we are considering. This does not necessarily 
mean that the ultimate theory of the electron, including quantum gravity, will attribute 
a timelike singularity to it, but at least this seems like an intriguing possibility. 

Another moral from our model is that timelike singularities are not as absurd as 
one might conclude from classical physics. Quantum mechanics (as well as Bohmian 
mechanics) differs in this respect from classical mechanics, which breaks down at timelike 
singularities because anything could come out of a timelike singularity, and there is no 
good reason for postulating any law about when and where particle world lines begin 
at the singularity. This is different in Bohmian mechanics because, as soon as we have 
defined the evolution of a wave function in Fock space, the probability current given in 
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terms of the wave function, such as 

J'^ = ?7''^, (1) 

defines the probabihty of particle creation at the singularity. Indeed, the stochastic law 
for particle creation is uniquely fixed if we strive for the following properties: reversibil- 
ity, equivariance, and the Markov property. Thus, the Bohmian particle positions are 
governed by two laws: an equation of motion, and a law prescribing the stochastic cre- 
ation of new particles. Both Bohmian and classical mechanics become indeterministic 
in the presence of timelike singularities, but in different ways: while Bohmian mechanics 
becomes stochastic, classical mechanics becomes lawless. One could say that our stochas- 
tic law closes exactly the gap in the laws caused by the timelike singularity. A timelike 
singularity causes concern in classical mechanics, but not in quantum mechanics. 

1.2 Horizons 

While the Reissner-Nordstrom example does not contain a horizon, we have a few 
remarks about the status of horizons in Bohmian mechanics. 

Consider a many-particles quantum system in a background space-time containing 
a black hole. From the point of view of orthodox quantum mechanics, it is natural to 
trace out all degrees of freedom inside the event horizon. From the Bohmian viewpoint, 
in contrast, this is not natural. Instead, it is natural to ask what happens behind the 
horizon. This is so because from the orthodox viewpoint, the most important question is 
what an observer will see, and it is a frequent assumption that the relevant observers sit 
at infinity. From the Bohmian viewpoint, the most important question is what actually 
happens. Thus, to define Bohmian mechanics in a curved space-time, we need to define 
as well the trajectories inside the black hole. (A viewpoint that, like the "Copenhagen" 
view of quantum mechanics, dismisses any theory of particle positions because it regards 
the latter as "hidden variables," may naturally tend to dismiss as unreal also everything 
that is hidden behind a horizon.) 

Thus, the Bohmian viewpoint leads us to the following attitude: What happens 
inside a black hole can and should be described by a physical theory. 

Indeed, by the nonlocal character of Bohmian mechanics, the velocities of the parti- 
cles outside the black hole may depend on the positions of the particles inside the black 
hole. But all this requires no additional research for the definition of the theory, as the 
Bohm-type law of motion we use (see Eq. (jl]) below) automatically implies influences 
across event horizons. What requires further work is not the presence of a horizon, but 
the presence of a singularity. (The need for this further work would disappear if none 
of the 3-surfaces belonging to the time foliation touched the singularity. However, for 
our purposes that is the uninteresting case. Furthermore, I see no reason why the time 
foliation should avoid the singularity.) 
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1.3 Concrete Model 



The concrete model we develop is based on the Dirac equation and the Reissner- 
Nordstrom geometry, but we expect that similar models can be constructed on the basis 
of other wave equations (such as the Weyl equation for massless right-handed spin-i 
particles, or equations for higher spins) and other space-times containing timelike singu- 
larities. We postpone the discussion of spacelike singularities to a future work; lightlike 
singularities behave, for our purposes, like spacelike ones. (A general singularity may 
well, like a general 3-surface, have different regions in which it is timelike, spacelike, and 
lightlike. Correspondingly, our separate studies of timelike and spacelike singularities 
apply separately to the timelike and the spacelike subset of the singularity.) 

The Dirac equation is known to have solutions of negative energy that are considered 
unphysical; indeed, the Hamiltonian is unbounded from below. We will not worry about 
this here, but rather treat the negative energy states as if they were physical. However, 
I recognize that, in a further development of this theory, the evolution should be so 
defined that only those states normally considered as physical can arise (in return, it 
may involve pair creation) . A difficulty one encounters when trying to write down such a 
model is the question: Which subspace exactly of the (one-particle) Hilbert space should 
be regarded as the space of physical states? Since the gap of 2mc^ in the spectrum of 
the Dirac Hamiltonian may disappear in the presence of an external gravitational (or, 
for that matter, electromagnetic) field, there is no longer any natural, gauge invariant, 
Lorentz invariant way of splitting Hilbert space into positive and negative energies. So 
we leave this question open. 

We also have to leave open the question whether the Hamiltonian we suggest here 
for the evolution of ip is truly self-adjoint; but there are reasons to believe this. If 
this proves right then our model will be an example of a theory of particle creation 
without an ultraviolet cut-off. In other words, in this setting the ultraviolet divergence 
which usually plagues any quantum theory with particle creation and annihilation will 
be absent, in fact without any special procedure a la renormalization. This encourages 
hope that the extreme curvature near a point mass may actually resolve, or at least 
ameliorate, the ultraviolet problem of QED and other quantum field theories. 

1.4 Motivation 

The investigations in this paper have several motivations, and lead to gains, in several 
areas: 

• It is a natural part of the research program on Bohmian mechanics to extend the 
theory to more general quantum theories, to all kinds of settings. To the extent 
that we have reason to believe that singularities exist in our universe, we obtain 
here a more appropriate version of Bohmian mechanics. Concerning techniques 
of constructing Bohm-type models, we find that the Bohm-type law of motion 
proposed by Diirr et al. [7] for relativistic space-time with a foliation works un- 
problematically even under the extreme conditions near a singularity; that not 
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every spacelike foliation is equally good for that purpose (as we discuss in Sec- 
tion [7]); that the stochastic approach to particle creation developed in "Bell-type 
quantum field theories" [TU] arises naturally in this context; that particle annihi- 
lation may well be deterministic and particle creation stochastic, even though the 
evolution is reversible. 

• Since Bohmian mechanics is a particularly precise and unambiguous version of 
quantum mechanics it may serve as a tool for studying quantum mechanics in 
curved space-time. Thus, our study can as well be regarded as one on quantum 
mechanics at space-time singularities. Concretely, we obtain a novel unitary evo- 
lution of the wave function in the presence of a timelike singularity, based on the 
Dirac equation and a boundary condition at the singularity (independently of the 
Bohmian viewpoint). 

• Several interesting features arise in the theory we develop: Timelike singularities 
emit particles; we can specify the law for how they do it; certain traits of quantum 
field theory (such as Fock space) arise from the presence of singularities; timelike 
singularities need not be unphysical. Finally, an ultraviolet regularization may not 
be necessary, thus suggesting a way out of the ultraviolet divergence for quantum 
field theories. 

I emphasize that the theory developed in this paper does not merely add Bohmian 
trajectories to otherwise known "orthodox" quantum theories. Instead, the Bohmian 
viewpoint helps define what the evolution equation should be, also for the wave func- 
tion. Indeed, the Fock space evolution corresponding to particle creation at a timelike 
singularity could be, but has not been, considered in an orthodox framework. The reason 
why it has not been considered before is, I think, twofold: First, since orthodox quan- 
tum physics tends not to ask what actually happens but rather what observers see, it is 
common to trace out the degrees of freedom inside a black hole, which means to ignore 
the phenomena this paper is concerned with. Second, in orthodox quantum physics it 
is common to focus on the scattering regime, i.e., on the asymptotic distribution of par- 
ticles a long time after their interaction. For defining a Bohm-like theory, in contrast, 
it is essential to define the trajectories, which may well depend on the particles hidden 
behind the event horizon, and this forces us to pay attention to what the time evolution 
laws are. Note also that a scattering theory can only be formulated in an asymptotically 
flat space-time, while our approach does not require asymptotic flatness or stationarity, 
even though the Reissner-Nordstrom geometry happens to be asymptotically fiat and 
stationary. 

The paper is organized as follows. In Section [2] we recall the relevant version of 
Bohm's law of motion. In Section [3] we recall the Reissner-Nordstrom geometry. In 
Section H] we discuss the general framework of particle creation and annihilation at a 
timelike singularity. In Section [5] we describe the concrete model. In Section [6] we 
draw parallels between our model and Bell-type quantum field theories. In Section [7] we 
describe a property of the time foliation that we have used. 
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2 Bohmian Mechanics 



2.1 Why Bohmian Mechanics? 

Bohmian mechanics is one of the most promising candidates for how quantum phe- 
nomena work. While the quantum formahsm describes what observers will see, and 
while the conventional "Copenhagen" view of quantum mechanics remains vague and 
paradoxical, Bohmian mechanics provides a possible explanation of the quantum for- 
malism in terms of objective events, in fact by postulating that particles have actual 
positions and hence trajectories. Bohmian mechanics is well understood in the realm of 
non-relativistic quantum mechanics, but needs further development in the directions of 
relativistic physics, quantum field theory, and quantum gravity. This paper concerns the 
relativistic extension in a classical gravitational field, but connects also with quantum 
field theory. 

The nature of quantum reality remains controversial. The traditional Copenhagen 
view is vague about what is real and what is not, and leaves the process of observation 
unanalyzed. To be sure, the predictions for the possible results of quantum experiments 
and their probabilities are uncontroversial. However, we are dissatisfied with the Copen- 
hagen view and feel that the quantum mechanical predictions call for an explanation. 
A good theory should provide an account of what objectively happens. Whether this 
account is true may be hard or impossible to decide, but we should require that it be 
consistent, in itself and with the quantum mechanical probabilities, while avoiding the 
vagueness of the Copenhagen version of quantum mechanics. 

Such an account is provided by Bohmian mechanics [H O [5], [11]; see [I3] for an 
overview. This theory postulates that particles have trajectories, governed by an equa- 
tion of motion of the type 

where Qt is the position of the particle at time t (or, for a system of several particles, 
the configuration), and j"^ and are, respectively, the quantum mechanical probability 
current and probability density as determined by the wave function ip. For example, if 
if) is a. Dirac wave function then 

f = ^^a^ , / = V'V • (3) 

As a consequence of the structure ([2]) of the law of motion, if at any time t the particle 
position (or configuration) is random with distribution p^* , then this is also true of any 
other time t. This property is called equivariance. As a (quite non-obvious) consequence 
of that, inhabitants of a Bohmian universe, consisting of these particles with trajectories, 
would observe the same probabilities in their experiments as predicted by the quantum 
formalism [IT]. That is how Bohmian mechanics explains quantum mechanics. In fact, 
Bohmian mechanics accounts for all phenomena of non-relativistic quantum mechanics. 

It is one of the virtues of the Bohmian viewpoint that it forces us to face those 
difficulties that can easily be swept under the rug in the orthodox framework. 
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2.2 Bohmian Mechanics and Relativity 

With the invocation of a preferred fohation T of space-time into spacehke 3-surfaces, 
given by a Lorentz invariant law and called the time foliation, it is known [5l [T71 [71 [25] 
that Bohmian mechanics possesses a natural generalization to relativistic space-time. 
The possibility of a preferred foliation seems against the spirit of relativity (see [17] 
for a discussion), but certainly worth exploring. It is suggested by the empirical fact of 
quantum non-locality, and it is suggested by the structure of the Bohmian law of motion 
([2]) for many particles, in which the velocity of a particle depends on the instantaneous 
position of the other particles. (But see [T3| for an indication of how a relativistic Bohm- 
like theory might be able to dispense with a preferred foliation.) Using a time foliation 
JF, a Bohm-type equation of motion was formulated in [7] for flat space-time (based 
on earlier work in |5j), and the straightforward generalization to curved space-time was 
formulated and mathematically studied in |25j : 

_^ cx r-^- (Xi(S), . . . , X^(S)) n n,^ (X,(S)) , (4) 

where Xfc(r) is the world line of particle k G {1, . . . , N}, t is any curve parameter, E 
is the 3-surface in JF containing Xfc(r), n{x) is the unit normal vector on E at x G E, 
Xj(E) is the point where the world line of particle i crosses E, and 

jW-Miv ^ ^(^^m (g, . . . (g) (5) 

is the probability current associated with the X-particle Dirac wave function ip. This 
wave function could either be a multi-time wave function defined on (space-time)^, or, 
since we never use ip for configurations that are not simultaneous, it suffices that ip is 
defined on the 3N + 1-dimensional manifold Usgjf '^^ simultaneous configurations. 

This generalization does not automatically include, however, space-time geometries 
with singularities. The treatment of singularities requires, as we point out in Section [H 
below, some fundamental extensions of Bohmian mechanics, and forms a test case for 
the robustness of the equation of motion (jlj). 

As mentioned before, the foliation might itself be dynamical. An example of a 
possible Lorentz invariant evolution law for the foliation is 

V^n^, - VuU^ = , (6) 

which is equivalent to saying that the infinitesimal timelike distance between two nearby 
3-surfaces from the foliation is constant along the 3-surface. This law allows to choose an 
initial spacelike 3-surface and then determines the foliation. A special foliation obeying 
([6]) is the one consisting of the surfaces of constant timelike distance from the Big Bang. 
Note, however, that the law of motion (jlj) does not require any particular choice of law 
for the foliation, except that the foliation does not depend on the particle configuration 
(while it may depend on the wave function). Note further that in a space-time with 
horizons, a foliation law like ([6]) will frequently lead to 3-surfaces lying partly outside 
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and partly inside the horizon. It will turn out from our study of singularities that ([6]) 
is not a good law for the time foliation in the presence of a timelike singularity as the 
foliations obeying (P) cannot be extended onto the singularity; see Section [7] for details. 



3 Reissner— Nordstrom Geometry 

Our main example of a timelike singularity will be the one of the Reissner-Nordstrom 
space-time geometry [191 ISSl [271 1201 [HI [IS], a solution of the coupled Einstein and 
Maxwell equations for a charged point mass. It is given by 

ds^ = -X(r)df + -^dr^ + r\dd^ + sin^ M^^) (7) 
A(r) 

with 

A r =1 + - 8 

and parameters M > and e G M. The Reissner-Nordstrom metric is spherically 
symmetric, static, asymptotically flat, and singular at r = 0. It arises as the gravita- 
tional field of a (non-rotating) point particle with mass M and charge e. We assume 
the so-called super- extremal case |e| > M; in this case A(r) > for all r > 0. In this 
context we note that if an electron were a classical point particle, its gravitational field 
would be a super-extremal Reissner-Nordstrom field containing a timelike singularity. 
This does not mean that electrons must, but does suggest they might, contain timelike 
singularities in nature. 

We will find it more useful to regard the singularity {r = 0} not as a line (as suggested 
by the picture of a point particle, and by the fact that in Minkowski space-time {r = 0} 
is a line) but as a surface (of topology R x §^). Thus, the space-time 

^ ^ {(t,r,cj) : t e M,r > 0,cj G = M X [0,oo) X (g) 

(where = means diffeomorphic) is a manifold with boundary. (For the mathematical 
notion of a manifold with boundary, see [161 [I] • We will sometimes write uo for a point 
on S^, and sometimes use the standard parameterization (■(9, ip) by polar and azimutal 
angle.) The singularity is the boundary 

^ = {r = 0} = (9^ ^ R X {0} X §^ (10) 

while the metric gfj,i, is defined on the interior = M x (0, oo) x S^. 

As the time foliation T for defining the Bohmian trajectories we use the foliation 
provided by the time coordinate t, i.e., the family of level surfaces of the t function: 

J^={E, :sGM} with = {x G ^ : t(x) = s} . (11) 

Concretely, = {s} x [0, oo) x S^. The fact that we are considering a static space-time 
invites us to regard S := [0, oo) x §^ as "the" 3-space. 
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4 Causal Trajectories at a Timelike Singularity 



Bohmian trajectories, the integral curves of the probabihty current, are causal curves, 
i.e., their tangent vectors are always timelike or lightlike (when based on the Dirac 
equation). Such a curve may well reach the singularity (indeed, in finite coordinate 
time t and finite proper time). What should happen if we consider a system of 
particles and one of them hits the singularity? The simplest possibility, or at least 
a natural possibility is that this particle gets annihilated, that it stops existing, and 
history proceeds with — 1 particle. This is the possibility we will study in this paper. 

In Bohmian theories, a configuration q of N identical particle is described by a set 
of N points in 3-space; thus, if S is 3-space, the configuration space is 

Q,v = {gCS:#g = iV}. (12) 

The wave function determining the motion of particles is best regarded as a function 
on Qn, which can be taken to be the covering space of Qn, or, somewhat simpler, 
Qn = S^. If the number of particles can change then we need several wave functions, 
one for each particle number. This leads us to considering state vectors in Fock space 

oo 

^ = 0^^, (13) 

N=0 

where J^fjsf is the Hilbert space of A^-particle wave functions, J^n = S±J^®^ , where 
^-i. means the symmetrizer (appropriate for bosons) and S*. the anti-symmetrizer (for 
fermions). Every element ip G can thus be regarded as a sequence 

^= (^o,V'i,---,V'iv,---)) (14) 

where at is an A^-particle wave function, defined on S^. Corresponding to the structure 
of Fock space, the configuration space is 

oo 

Q = U 2iv, (15) 

where Qat is the configuration space of particles, as defined in (fT2|) . The Fock state 
ip can then be regarded as a function on (the covering space) 

oo 

Q=[jQN. (16) 

N=Q 

(Both Qo and Qq contain just one element, the empty configuration.) 

Thus, in the theory we are developing, if one among A^ particles hits the singularity, 
the configuration jumps from Qn to Qn-i in a deterministic way: q f{q), where the 
jump function / : OQn U„<ArQ„ is given by 

/(g) = g\9S. (17) 
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Here, dH, the boundary of S, is the intersection of E with the singularity S = d^; 
OQn is the set of configurations in which at least one particle is located on 9E; and 
ffT7|) means that / just removes all particles located on 9E, i.e., all particles that have 
arrived at the singularity. 

The situation in which the configuration space is a manifold with boundaries has been 
considered before in [12]; we recapitulate the Bohmian dynamics of the configuration 
developed there. It includes deterministic jumps that occur whenever the configuration 
reaches the boundary dQ of configuration space. We desire three properties of the 
theory: reversibility, equivariance, and the Markov property; as pointed out in |12j . 
these properties entail stochastic jumps with a certain rate. 

This leads us to a stochastic process {Qt)tm. in Q consisting of smooth motion 
interrupted by jumps from the boundary to the interior or vice versa. The jumps from 
dQ to Q° are deterministic and occur whenever the process hits the boundary. The 
jumps from Q° to dQ are stochastic, and their rates are fully determined by requiring 
that (i) the process is Markovian and equivariant, and (ii) the construction is invariant 
under time reversal, in that the processes associated with suitably time-reversed wave 
functions are reverse to each other, in distribution. 

The deterministic jump law, abstractly 

= f{Qr-) (18) 

for a fixed mapping / : dQ Q°, is in our case given by (1171) . Since we want the 
theory to be reversible, we must also allow for spontaneous jumps from interior points 
to boundary points. Since we want the process to be an equivariant Markov process, 
the rate for a jump from q' G Q° to a surface element dq (1 dQ must be, as one can 
derive, 

'yt{Q'^dq) = ^^^^uidq,q'), (19) 

where = max(s, 0) denotes the positive part of s G M, is the probability current 
relative to coordinate volume (as discussed in detail in Section 15.11 below), is its 
time component and thus the probability density relative to coordinate volume, f is its 
component orthogonal to dQ (in our case, the radial component), and v{dq,q') is the 
measure- valued function defined in terms of the volume measure n on Q and the surface 
area measure A on (9 Q by 

for every B C dQ. In our case, the measure is on Qm locally the product of copies 
of the coordinate volume measure on E = [0, cxd) x S^, and the measure A is on dQ^ 
locally the product of the surface area measure on §^ and — 1 copies of the coordinate 
volume measure on E = [0, oo) x S^. Note that the set of configurations with two or 
more particles on the singularity has A-measure zero. It follows that z/(-, q) is a measure 
concentrated on the set of those configurations that consist of q plus one further particle 
on the singularity (9E = S^, and the measure is essentially a copy of the surface area 
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measure on S^. Thus, concretely, the creation rate density relative to the area measure 
on §^ is 

at{q-^qUu) = — — , 21 

where we have simply written u for the point (r = 0, u) on the singularity S" fl S. 

Note that a jump to a boundary point q at which the current is pointing towards 
the boundary, f{q,t) < 0, would not allow any continuation of the process since there 
is no trajectory starting from q. The problem is absent if the velocity at q is pointing 
away from the boundary, f{q, t) > 0. (We are leaving out the case _f (g, t) = 0.) On the 
other hand, jumps from q to f{q) cannot occur when j{q,t) is pointing away from the 
boundary since in that case there is no trajectory arriving at q. Thus, the jumps must 
be such that at each time t, one of the transitions q — > /(g) or /(g) — g is forbidden, 
and the decision is made by the sign of f{q, t). 

Given the law of motion 

dt f 

together with the deterministic jump law fllSp and stochastic jumps with rate (1191) . we 
obtain the following probability transport equation at g' G Q°: 



^{q',t) = -d,f{q\t) - at{dQ,q')f{q\t) + j u{dq,q'] 



liq.ty. (23) 



dQ 



For equivariance we need that (l23i) . has the structure of the transport equation for J* 
that follows from the Schrodinger equation. This is what our model achieves by means 
of a suitably modified Dirac evolution and a boundary condition at the singularity. 



5 Concrete Model 
5.1 Probability Current 

When we speak of probability density, or probability current density, it plays a role that 
densities can be considered either relative to coordinate volume or relative to invariant 
volume (proper volume), and this leads us to considering two current vector fields on ^ 
(and afterwards on Q), which we denote respectively and j^. We elucidate this below 
in more detail, beginning with the single-particle case. The reason we consider densities 
relative to coordinate volume is that they stay bounded at the singularity, whereas the 
relevant densities relative to invariant volume diverge. 

The current given by a one-particle wave function, = ip'^^^ip, is the current relative 
to invariant volume, which means the following. For every spacelike hypersurface E C 
let n'^(x, E) be the future-pointing unit normal vector on E at x G S; let \Y.{d^x) 
be the 3-volume measure associated with the Riemannian 3-metric on E inherited from 
g^^j. The Bohmian trajectory associated with ^ is a randomly chosen integral curve L 
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of whose probability distribution P is such that, for every spacehke hypersurface S, 
the trajectory's unique intersection point with E has distribution j^n^Xj^, i.e., 



P(L n S C A) = / fix) g^,{x) rf^x, S) As(c/'x) 

J A 



(24) 



where A is any subset of S, L n S the (set containing only the) intersection point, and 
P(L n S C A) the probability of this point lying in A. 

When computing densities relative to coordinate volume, we regard the coordinate 
space, in our example M x [0, oo) x as a Riemannian 4-manifold ^ with metric g^j^u, 
in our example 

Every spacelike hypersurface S C ^ can also be regarded as a hypersurface in let 
ri'^(x, S) be the vector at x G S that is normal in coordinates (i.e., in g^^,) to S, has 
unit length in coordinates, and points to the same side of S as n^{x, S). Likewise, let 
\Y,{d^x) be the 3- volume measure defined by the 3- metric inherited from g^^,. As we 
show presently, there is a vector field such that, for every spacelike hypersurface S, 

r{x) ~g^, n\x, S) \^{Sx) = j^(x) g^^ix) n^(x, S) Xj,{d'x) . (26) 

Indeed, 

r(x) = ^^ix)fix), (27) 

where e((i^x) is the invariant 4- volume measure associated with g^^, and e{d^x) the co- 
ordinate 4- volume measure (associated with gf^^)- Note that, as a consequence of the 
proportionality, j'^ and j'^ have the same integral curves (up to an irrelevant reparame- 
terization) . 

To show (!26|) and fl27j) . we use differential forms (see, e.g., [H]). Let j be the 3-form 
obtained by inserting into e, the invariant volume 4-form associated with g^^,. Then 
the distribution (EH) can be written as 



P(L n S C A) = [ j 

J A 



(28) 



for every subset A of S. That is, j encodes the probability current without reference to 
any metric. Making the same step backwards with g^^,, we obtain that 

^fiunX 3 vkX ^fiUKX ; (29) 

which implies f l27p . 

One can compute the Radon-Nikodym derivative between the two measures by 

(k ^ y- det c/^^ ^^^^ 
de ^- det g^,^ 
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In our example case, de/de = r^. (Note that taking (0, oo) x as coordinate space 
differs by a factor of from taking \ {0}.) We set 




de 
de 



(31) 



so that = t/j^'^ip. 

If E is not spacelike, the distribution of the intersection point L n E is still given by 



for any ^4 C E, provided that L intersects E at most once, and that always points 
to the same side of E as j^. In particular, if E = S* is the singularity then we obtain 
that the probability distribution of the random point where L hits the singularity has 
density (in coordinates) equal to (the absolute value of) f, the radial component of J'*, 
except when _f > so that there is a current away from the singularity. 

Since we only need that j'^ remains bounded at the singularity, we can allow j'^ to 
diverge like r~^. This means, since = i/^^'^i^, that ip can diverge like while 7^* 
stays bounded, or 7^ can diverge like while ip stays bounded. We choose the second 
option, so that relevant wave functions ip remain bounded at the singularity. This is 
connected to the definition of the spin spaces at the singularity. 

5.2 Spin Spaces at the Singularity 

The spin spaces '^^^ standardly defined at points x G , and we also need to define 
them at points x E S in order to be able to talk of -0(0;) and j^{x) for x e S. As a 
preparation, let us look at how the tangent spaces are defined for x E S. (We find 
it useful to regard as a vector space rather than half-space.) Using coordinates, 

we can simply say, for x — {t,r — 0,u!), that 



where the first factor M is understood as the tangent space to the manifold M at t, and 
the second as the tangent space to the manifold IR at 0. 

As a part of this specification of tangent spaces, we intend a certain topology (and 
different iable structure) on the bundle = Uxe^T^^: for example, if we set x = 
(t, r,uj),y = (t, 0, u), keep t and uj E E>'^ fixed and let r — > so that x y, then the unit 
normal to the level surface E^ of the t function, n^{x), converges to 0, as its coordinate 
components arc (A(r)~^/^, 0, 0), while the unit normal relative to g^j,^, h^{x), converges 
to h^{y). Thus, another way of defining the tangent spaces on the singularity consists 
in specifying a tetrad field in (i.e., a basis for every tangent space) that possesses 
a smooth extension to d^. In our example, we can take as the tetrads in simply 
the coordinate bases (or, orthonormal bases relative to g^u)- 

In the same way, we define the spin space for x e by specifying a ba- 
sis bx of yx for every x e and postulating that the field bx possesses a smooth 



P(L n E C A) = / r{x) ~g^,,{x) h%x, E) As(rf'a;) 



J A 



(32) 



Tx^ = R X M X T^S^ , 



(33) 
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extension to . To this end, fix any x G and begin with the coordinate ba- 
sis e-x of Ta;^, = {.dt,dr,d^,d^p). By normahzing these four vectors, we trans- 
form them into an orthonormal basis (i.e., Lorentz frame, a.k.a. Minkowski tetrad) 
= (A^^^^^t, \^^'^dr,r^^d^, {r sm'd)~^dip). To this orthonormal basis there corresponds 
a basi^ bx of S^x] the correspondence is canonical up to an overall sign which we 
choose continuously in x; is an orthonormal basis relative to the scalar product 
07^ (x) g^u{x) n^i^x) ip. This basis we rescale by \{r)~^/'^ [dej de)~^l'^ to obtain b^- This 
completes the definition of b^- 

In other words, if we define the sesquilinear mappings a, a : x 5^x — > CT^^ by 

a(0,^) = 07'^^, a(0,^) = 07'^V>, (34) 

then the coefficients of a relative to bx and tx are the usual three Dirac a matrices 
together with oP = I, the identity matrix, while the coefficients of a relative to bx and 
follows: 

{a\ a\ a\ = (^/, a\ ^a', -^^^j^^—^a'^ , (35) 

where the a* on the right hand side denote the standard Dirac a matrices. This equation 
could also be taken as the definition of bx- 

From fl35|) it follows that the coordinate formula for the probability current is given 

by 

/ = ^*a''^ (36) 

where the components of are taken relative to the basis bx, and the matrices are 
given by the right hand side of (l35l) . Note that these matrices possess the following 
limits 

asr^O: (a*, a^ a'', a^) = fo, a\ -a^, -a^V (37) 

V e esmty / 

with e the charge parameter in (J?!). As a consequence, the probability current at the 
singularity always has vanishing time component, i.e., is tangent to S^. 



5.3 Boundary Condition 

We impose a boundary condition at the singularity on the wave function. Let 2; be a 
point on 5* n S and xi, . . . , xn-i G S \ 5. A Dirichlet boundary condition [23] 

i){xi, . . .,xn-i, z) =Q (38) 

would lead to the (for our purposes) uninteresting behavior that no particle ever reaches 
the singularity. The boundary condition that we want should allow that there can 
be a probability current into the singularity but make sure that that current into the 
boundary of Qat corresponds to the gain in probability in Qn-i- 

-'^The same can be done in two-spinor calculus, where this basis is called a "spin frame" |21j . 
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To this end, we choose two functions (f)± : Ei^ —>■ such that for every G S^, 
4>±{u>) is a normahzed eigenvector of with eigenvalue ±1 (recall that the Dirac a 
matrices are Hermitian 4x4 matrices with eigenvalues ±1). Being eigenvectors to 
different eigenvalues, the two are orthogonal, 0+(c<j)* 0„(co') = 0. As a consequence, 0± 
are two orthogonal vectors in L^(§^,C^), each with norm a/Itt. Thus, the 0± span a 
2-dimensional subspace W C L^(S^,C^). In the following we will simply regard the 
singularity S' fl S as and write u for a point on the singularity. 

We present the boundary condition as the conjunction of two boundary conditions, 
the first of which reads 

for every q E Qn- Equivalently, il){q,uj) can be written as 

iP{q, u) = (^c+{q) M^) + c^q) </)_(aj)) ® xiq) (40) 

with complex coefficients c±{q) and suitable x(g) G (C'')®^. The second boundary 
condition reads 

V^(g) = ^(c+(g) + c_(g))x(g), (41) 
Voir 

where x{q) denotes the same spinor as in fHOj) . Equivalently, we can write fHT]) as 

^{q) = ^{Moo) + <j>4u;)y^{q,u;) G §^ (42) 

VoTT 

So we have specified what the boundary conditions are. They are obviously linear. 

Technically speaking, the boundary conditions are part of the specification of the 
Hamiltonian H, as to specify H means to specify (i) its domain (a dense subspace of 
Hilbert space J^) and (ii) how H acts on vectors in the domain. The domain consists 
of those wave functions from the first Sobolev space (i.e., possessing a square-integrable 
first derivative in the sense of distributions) that fulfill the boundary conditions described 
above. On such a function H acts as follows: For g G Q°, 

Hi/jiq) = -iha\q)-^a\q) di^{q) + V{q) ^{q) + Hj ij{q) , (43) 

where the index i runs through 3A^ dimensions if g G Qn, and V{q) assumes values 
in the Hermitian matrices on (C^)®^ and includes the mass term (3m, all connection 
coefficients that may arise from covariant derivatives, the external electromagnetic field 
(if desired), and all other potentials. The first two terms are just the (appropriate version 
of the) Dirac equation, while the third term creates a link between ip^ and ipN+i'- 

Hiip{q) = ihV2^ I {(p+iuj) - 0_(cu)) V(g,cu) du = (44) 
= thV32^{cM - c_(g)) x{q) ■ (45) 



(39) 



15 



What needs to be checked is the conservation of probabihty. At {g U : G S^}, 
the probabihty current into the singularity is 



f{qlJuj)duj= / 'ijj{q,ujy a'^{q,u) ■ijj{q,uj) duj 



(46) 



^|J{q, u)* a;*(g) ® ip{q, uj) doj 



(47) 



n A(r(x)) 

xiqYxiq) 
n A(r(x)) is 



c+{q) + c^{q) 0-(^) ) tt^ ( c+(g) 0+(^) + c_(g) 0-(^) ) c/w 



(48) 

c+(g)0+(u;) + c_(g)0_(a;)) fc+(g) 0+(u;) -c_(g) 0_(u;)) rfu; = (49) 



47r 



X(9)*X(9) 
n A(r(a;)) 



c+(?)r-ic-(?)r)- 



(50) 



On the other hand, since the Dirac equation imphes the conservation of probabihty, the 
gain of probabihty at g G Q° is the contribution to the increase of j*(g) due to Hi: 



dt 



Hi 



ilm^{qya'Hj^/j{q) 



A(r(x)) 

x£q 



Imiplq)* Hjiplq) 



(51) 



4vr||^^ Re K(g) + c_(g))*K(g) - c_(g)) 

xGq 

n A(r(x)) V 



(52) 



(53) 



which is the same as (1501) . This shows that probabihty is conserved. We thus beheve 
that (l43l) . together with the boundary conditions, defines a self-adjoint Hamiltonian and 
thus a unitary evolution on Fock space. 



6 Comparison to Bell- Type Quantum Field Theo- 
ries 

Two ways of extending Bohmian mechanics to quantum field theory are known: either by 
postulating that a field configuration (rather than a particle configuration) is guided by 
a wave function (understood as a functional on the field configuration space) [HE], [23], or 
by introducing particle creation and annihilation into Bohmian mechanics [2], [9], [TOl [12] 
(but see also [6] for a third proposal). The second approach is called "Bell-type quantum 



16 



field theory," as the first model of this kind (on a lattice) was proposed by Bell p]. 
In these theories, the motion of the configuration along deterministic trajectories is 
interrupted by stochastic jumps, usually corresponding to the creation or annihilation 
of particles. The jumps are governed by the following law prescribing the jump rate 
a"^ (probability per time) in terms of the wave function ip, which is usually from Fock 
space: 

, ._, [ir^m)mM')m)v ...^ 
{mm) ' ^ ^ 

where Hj is the interaction Hamiltonian. The configuration Qt thus follows a Markov 



process in the configuration space Q of a variable number of particles, as defined in (fTS!) 
above; see [10] for a detailed discussion of this process. 

We expect that further analysis will show that our model for stochastic particle 
creation at a timelike singularity fits into the scheme of Bell-type quantum field theories 
when understood in the appropriate way. In other words, we expect that the jump rate 
( PTj) can be regarded as a special case of 



7 Requirements on the Time Foliation 

In the presence of singularities, difficulties arise about the choice of the time foliation 
that had not been considered previously. While the foliation we used on the Reissner- 
Nordstrom space-time, given by the t coordinate, worked well for our purposes, there also 
exist problematic foliations. Here is an example: Start with the {t = 0} hypersurface 
and propagate it according to (E]), i.e., push every point of the hypersurface to the future 
at the same rate (so that the infinitesimal slice between two nearby hypersurfaces has 
constant proper thickness). The hypersurfaces we thus obtain are the level surfaces of 
the function T that yields, for a space-time point x, the timelike distance of x from the 
hypersurface {t = 0}. At first sight, this might seem like a reasonable foliation, but in 
fact it does not cover all of not even That is because the T function is infinite 
in a large region of space-time, in fact (in the future of {t = 0}) at all points on or later 
than the outgoing radial null geodesies starting at t = 0. Let me explain. 

I begin with describing the outgoing radial null geodesies of the Reissner-Nordstrom 
geometry: They are the integral curves of the vector field = \{r)~^/'^dt + A(r)^/^(9,., 
as this vector is null (easy to check), radial {u^ = = m*^), future-pointing (u* > 0), and 
outward-pointing {if > 0). Thus, the outgoing radial null geodesies are explicitly given 

by 

with constant d and (p. For small r, t{r) ^ r^/3e^. Correspondingly, the incoming radial 
null geodesies are given by 

r dr' 

'''■'^"'"-/o '''' 
Note that at every point of the singularity, one incoming radial null geodesic ends and 
one outgoing radial null geodesic begins. 
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Consider the curve fl55l) for any t{0) > and call it C. Any point a; on C must have 
T = oo: If it had finite T value, then choose ri so small that A(ri) > T{x)/t{0) and 
ri < r(x); that is possible because A(r) — > oo as r — > 0. Let Xi = {t{ri),ri,'d{x),(p{x)), 
which lies on C, too. Then the parallel to the t axis (in coordinate space) through xi is 
a timelike curve connecting xi to the hypersurface {t = 0}, and the proper time along 
this curve between xi and {t = 0} is A(ri)t(ri) > A(ri)t(0) > T(x). Since T{xi) is the 
supremum of such curve lengths, we have that T(xi) > T(x), so that x must lie in the 
past of the spacelike hypersurface {T = T(xi)}, in contradiction to the fact that x lies 
on the future light cone of xi. 

Since t(0) > was arbitrary, it follows that any point in the future of the hypersurface 
defined by fl55|) for t(0) = and arbitrary i), (p, must have T = oo. And thus, T foliates 
only a part of Moreover, the leaves cannot be extended in a disjoint way to the 
singularity, as any {T = const.} hypersurface has {(t = 0,r = 0)} x §^ as its boundary 
on the singularity. One could say that the leaves intersect on the singularity. 

In contrast, the t function provides a foliation of all of including the singular 
boundary. That is a property we need of the time foliation. 
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